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1. Introduction 

Bifurcation involving chaotic motions is one of the origins of diversity in complex 
systems and reconnection among stable/unstable manifolds is such an example. One 
of the simplest systems exhibiting reconnection is the Harper map pflEj- ^ or practical 
applications, it is helpful to have an analytical view and we analytically study the 
reconnection in the nearly integrable Harper map. 

To study weakly perturbed stable/unstable manifolds, the conventional Melnikov 
perturbation method is not applicable since splitting between stable and unstable 
manifolds is exponentially small for small perturbation parameter a. The method of 
asymptotics beyond all orders (ABAO method) is one of the useful methods dealing with 
such situations. The key idea is to employ the so-called inner equation, which magnifies 
the behavior of the solution near its singularities, and to apply the Borel transformation 
for investigating the inner equation and analytically continuing its solution. The 
idea of the inner equation was first used by Lazutkin, Schachmannski, TabanovjSj to 
derive the first crossing angle between the stable and unstable manifolds of Chirikov's 
standard map and by Kruskal and Segurjl] to study a singular perturbation problem of 
ordinary differential equations. Hakim and MallickjH] introduced the technique of Borel 
transformation within in this context and obtained the first crossing angle previously 
obtained in [Oj. Tovbis, Tsuchiya and Jaffe 7] improved their method and derived 
analytical approximations of perturbed stable/unstable manifolds for the Henon map. 
In these approaches, only the dominant part of the inner solution is taken into account. 
Later, Nakamura and Kushibe ^H] proposed a method of systematically improving 
the solution of the inner equation with the aid of the Stokes multiplier and studied 
Chirikov's standard map. The ABAO method was applied to some higher dimensional 
systems as welljTTJ [TJ]. Also, Gelfreich and his collaborators EH EEH1 EH EI are 
studying the splitting of separatrices and related bifurcations for various systems by a 
slightly different but more rigorous way. 

Roughly speaking, the procedure [7| l8| l9| ITf)] obtaining an approximate expression 
of the unstable manifold is summarized as follows: 

(i) Find singularities in the complex time domain of the lowest order ordinary 
perturbative solution of the unstable manifold. 

(ii) Magnify the neighborhood of a singularity closest to the real axis and derive 
an asymptotic expansion of the lowest order solution which is valid in a sector 
containing t = — oo. 

(iii) The method of the Borel transformation is used to construct the asymptotic 
expansion which is valid in the other sector. Usually, there appear additional terms 
which are exponentially small for real time values. 

(iv) Going back to the original equation, derive equations for corrections corresponding 
to the exponentially small terms found in the previous step. An appropriate solution 
is chosen by matching its asymptotic expansion with that obtained in the previous 
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step. 

The additional terms are determined by the singularities of the Borel-transformed 
solutions and, so far, only contributions from poles were considered. In this paper, with 
the aid of ABAO method, the reconnection of the unstable manifold is studied for the 
Harper map. Our analysis mainly follows the method by Nakamura and Kushibe |10| . 
but contributions from the branch cuts of the Borel-transformed solution are also taken 
into account. 

The Harper map depends on a real parameter k and is defined on (v, u) G [— 7T, n] 2 : 
v(t + a) — v(t) = — crsinti(t) 

u{t + a) — u(t) = ka sin v(t + a) (1) 

where a(> 0) is the time step and plays a role of the small parameter. Since the case 
of k < is conjugate to that of k > 0, it is sufficient to consider the latter In the 
continuous time limit, where a — > 0, the map reduces to a set of integrable differential 
equations: 

v'{t) = — sin u(t) 

u'{t) = ksinv(t) . (2) 

Throughout this paper, The prime is used to indicate the differentiation with respect 

to time t such as u'{t) = — ■• Eq.(0) admits topologically different separatrices 
depending on the parameter k (cf. Fig. and the separatrix changes its shape when 
k — > 1. Since the solution for k > 1 can be obtained from that for k < 1 by a simple 
symmetry argument (cf. |Appendix A| ), we restrict ourselves to the case of < k < 1 
and consider the change of the stable/unstable manifolds for k — > 1. In the text, 
an approximation of the unstable manifold is constructed and only the result of the 
approximate stable manifold is given. 

The rest of this paper is arranged as follows: In the next section, the lowest order 
solution of the Melnikov perturbation method is constructed and the inner equation is 
derived. In Sec. 3, the Borel resummation method is used to solve the inner equation and 




Figure 1. Separatrix of the Harper map in the continuous limit 
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derive the additional exponentially small terms. In Sec. 4, the solution of the original 
equation with the corresponding corrections is obtained and compared with numerical 
calculations. The reconnection of the unstable manifolds is discussed in Sec. 5. Sec. 6 
is devoted to the summary. 



2. Melnikov Perturbation and Inner Equation 

Let Vq,Uq be the solutions obtained by Melnikov perturbation, i.e., 



v 



(t) EE J2 ^ V 0n(t) , Uo(t) = J2 ^ U 0n{t) 



n=0 



n=0 



The lowest order terms satisfy (J2J) and the first and second order terms obey 

1 



v' 01 (t) + U l(t) COS U 00 (t) 



u' m (t) - kv m {t) cos v 00 it) = 2 U od(*) 



(3) 



11 „, u 2 (t) 
2 o 2 

u' 02 {t) - kv 02 (t) cos w 00 (t) = ttm'oi (*) - c m m(*) - ^-^TT— sin u oo(0 (4) 

2 2 

The lowest order solution of the unstable manifold with vq(0) = is given by 
v oo(t) = — 2 tan -1 



ttQo(i) = 2 tan 1 



\J\ — k sinh Vct 

~~ k~~ 1 

1 - k cosh v^t 



(5) 



Because of the boundary conditions v 0n (t) — > 0, Mon(^) - >■ for t — > — oo and t>o n (0) = 0, 
the first and the second order solutions are 



u oi(t) = -yi{t) 



v 01 (t) = , 

v 02 (t) = - ^ + ari(*) I ^ - 2(1 - fc) 

Mo2 (t) = 7T 1 7 2y[(t)- yi (t) I kt - 2(1 -k) 



(6) 



v^fc tanh Vkt 
1 — k tanh 2 Vkt 

V~k tanh Vkt 



1 — A; tanh 2 Vkt 



(7) 



where auxiliary functions are given by 
x x (t) = -2y/k(l-k) 



cosh Vkt 



y x (t) = -2ky/l-k- 



1 + (1 - k) sinh 2 Vkt 

sinh A/fct 
f (1 - k) sinh 2 Vkt 



(8) 
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Figure 2. Singular points (crosses) of the lowest order solution vqo, uqo in the complex 
time domain. 



We observe that the lowest order solution (jSj) has two sequences of singular points 
in the complex t-plane (see Fig. (2j). 



1 f l±Vk ( 1\ 



/ — I In ("■■-) ni ? (0) 

Now we examine the behavior of perturbative solutions near two singular points 
in the upper half plane closest to the real axis: t\ = |ln ^= + |7rz|, t<i = 

^= |ln + §7ri |. From now on, t c stands for ti or t 2 , and the upper sign corresponds 
to ti and the lower one corresponds to t 2 - 

As easily seen, the perturbative solution admits an expansion like 

00 00 (n) 

Mt,a)-v Q0 (t)=^^-^— l a" (10) 

n=l 1=0 V lc > 

This indicates that the higher order terms with respect to a has higher order poles at 
t = t c and that the behavior near the singularities should be carefully investigated. For 
this purpose, it is convenient to introduce a rescaled variable z = El El El ITU] . 
Then, one finds that & (z,a) = v (t c + az,a) — v 00 (t c + az) and ^ (z,a) = u (t c + 
az, a) — u 00 (t c + az) can be expanded into power series with respect to a: 



00 



1=0 1=0 
where ($ o,^oo), (^oij^oi) ••• correspond to the most divergent terms, the next 
divergent terms .... 

We remark that, in terms of §(z, a) = v(t c + az, a) — voo(t c + az) and ^(z, a) = 
u(t c + az, a) — u 00 (t c + az), the Harper map (JTJ) reads as 

A&(z,a) = — a sin \l> [z, a) cos u 00 (t c + az) 

— a cos \1/ [z, a) sin u 00 (t c + az) — Av 00 (t c + az) 
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A^f(z -l,cr) = ka sin $(0, a) cosf 00 (t c + az) (12) 

+ ka cos a) sinf o(t c + — A-u o(t c + &(z ~ 1)) 

where A stands for the difference operator: A/(z) = /(z + 1) — /(z), and that 
$ (z, cr), ^ (z,a) are its perturbative solutions with respect to a. This equation will 
be referred to as an inner equation. 

Substituting (fTT|) into (fT2"|) . one can easily obtain the equations for <3>oz, ^oz- The 
first two sets are as follows 



A$ 00 (z) 
A^ o(z) 
A$„i(«) 
AtfoiCz) 



e q=i*oo(z) / J 

z zln 1 + - 

z V z 



2 + 1 



± i In [ 1 + - 

z 



± 



1 | ^oi(z) 



3 =F«*oo(2) 



+ -(* 



$01(3 + 1) 



.1 - fc' 



2+1 

By matching the residues at z = 0: 

'$0fc(z) 
*0*(*)/ _ (J 



e i*oo(^+l) + ^ _ fc ) 



Res 



Res 



V k+l(t) 
U k+l(t) 



t=t c 



one obtains the solution corresponding to the expansion (|10j) 



^oo(^) 
^oo(^) 
$01(2) 
*oi(«) 



12z 2 



107i 
432(k 4 



O A 



T 7; 



+ 



191i \ 

2z 2Az 2 24z 3 ' 8640?/ 
% ±kt c + 1 
24 z 



+ 4 



+ 4 



i(k - 1) i k(±t c + 1) i k(±t c + 1) 



+ 4 



(13) 



(14) 



(15) 



(16) 



;i?) 



4 24 z 48 z 2 

This is the asymptotic solution of the inner equation which is valid in the sector 
involving Re z = — 00. It is necessary to analytically continue it to another sector 
involving Re z = +00 in order to construct an approximation of the unstable manifold 
in the time domain Re t > Re t r . This will be carried out in the next section. 



3. Borel Transformation and Solution of Inner Equation 

3.1. Borel Transformation and Analytic Continuation 

As mentioned in the previous section, in order to construct an approximation of the 
unstable manifold in the time domain Re t > Re t c , the solution of the inner equation 
in the sector involving Re [z] = —00 should be analytically continued to the sector 
involving Re[z] = +00. The analytic continuation is carried out with the aid of the 
Borel transformation. Let us begin with Q 00 (z), ^qq(z). 
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Figure 3. Integral path in p-planc 



Since we are interested in the unstable manifold, we define the Borel transforms 
Vo(p), U Q (p) of $00(2), ^00(2), respectively, as 



p—oo 

$oo(z) = L[V (p)](z) = / dp e~ pz V (p) 

Jo 

* 00 (z) = L[U (p)](z) = / dp e~ pz U (p) 

Jo 



The analytic continuations of $oo(z) and ^oo(z) from Re[z\ < 0, lm[z] < to Re [2;] > 
are obtained by rotating the integral path in the p-plane from C to C (cf. Fig. |BJ). 
Other possible singularities and branch cuts of Vo, Uq are also shown in Fig. EJ Note 
that the integral path should be rotated over | [rad] in a clockwise way so that the 
convergence domain rotates counterclockwise. Thus, the analytic continuation from 
Re[z] < 0, lm[z] < to n + e < argfz] < 2ir + e is given by 

/ dp e~ pz V (p) -> / dp e~ pz V (p) - / dp e~ pz V (p) (19) 

JO Jo J-y 

where e is a small positive real number. 

If one can write down the Borel-transformed inner solution explicitly, one has only 
to change the integral path. But it is difficult and we have to deal with the inner 
equation and the Borel-transformed inner equation simultaneously. The procedure can 
be summarized as follows 

(a) Write down the equation for V , U : 

[ T i)-Ut\p) 1-e" 



/p _ _ 



n=l 



±i(l-e')t/„( P ) =1+ /" dA^^t^-^-K (20) 



" I n=l y 



where Vq*^ denote the nth convolution defined by 

V [ * n) (p)= [ P dxV (p-x)V { * {n - 1)] (x) (21) 
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This equation indicates that the Borel transforms Vq(p), Uq(p) may have 
singularities at p — ±2m and branch cuts starting from them (cf. Fig. EJ). This 
observation and the nonlinearity of the inner equation imply that additional terms 
like — 27rzRes (V (p)e~ pz ) \ p=27rin and those from branch cuts may appear in Re[z] > 0. 
(b) Observation (a) suggests the expansions 



— 2irinz 



n=0 

oo 

*(^)Uo = E $ «°W e ~ 2 ™ 2 ( 22 ) 

ra=0 

where <J> n o, ^ n o consist of finite polynomials of z and power series of 1/z. Solve 
the equations for $ n0 , ^ n o, which can be obtained by substituting (j22J) to (|T2|) and 
comparing the same order terms with respect to e~ 2mz . Note that $00; ^oo are the 
solutions of (fIH| . 



(c) With the aid of the formula 

- 2iri(-z) j ~ l e- 2ninz = - I dp e~ pz 
-2m 



■_ e -2ninz = _ / dpe ~pz 



7 



U - 1)' 

(p — limiy 
(p — 27rinY ln(p — 2irin) 



z?- 

find V (p), U (p) in the upper half p-plane from the following equation 



p OO p oo 

/ dp e~ pz V (p) = ®no(z)e~ 2 ™ nz , - dp e~ pz U (p) = ]T * n0 {z)e 

h n=l T n=l 



2-Kinz 



In the above, we choose the branch cut of the logarithm along the positive imaginary 
axis as shown in Fig. EH Note that this choice fixes the Stokes line on the negative 
imaginary axis in the z-plane, which corresponds to the line in the original t-plane 
joining t c and its complex conjugate. 

(d) By the procedures explained so far, only the singularities in the upper half p-plane 
are taken into account. The contributions from the singularities in the lower half 
plane are determined from the fact that Vo, Uq is pure imaginary when z is real. 

(e) Determine Vo, f/o by comparing the MacLaurin expansion of the analytical 
guess obtained from (a)-(d) and the numerical power series solution of the Borel 
transformed inner equation ()20|) . 

3.2. Determination of Vo, Uq 

In this subsection, following the steps (b)-(e) explained before, Vo, U are calculated. 
The equations of $io, ^10 are given by 

e =Fi*oo(«) e «*oo(2+l) 

A$ 10 (z) = ± *io(*) , Att 10 (*) = ±$ 10 (* + 1) — (23) 

Z -2 + 1 
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which have two independent solutions: 



(A) 
10 

(B) 
10 



1 23 

z + 48z J 



z + 



+ 



17 



(.4) 



10 



±|i 1 



1 



12z 1080z 3 



+ 



(B) 
10 



1 1 

± ' r+ 2 " 540^ 



+ 



and $io, ^io are their linear combinations. 



$10 
#10 



ca 



^10 
^10 



+ c B 



(B) 



From 



Vo(p) is given by 

„ CO 

- / dpe-**V (p)= J^^oe 

J T 1 



10 
^(B) 



2irinz 



(24) 



n=l 



cb\z + 



1 



-27TJ2 



(25) 



12^ J z \z" 

and Uq can be represented in a similar way. As easily seen from (|20|). the coefficients 
and cb are imaginary number. 

In order to determine the Borel transforms Vq, Uq from the expansion (|2*3jl by taking 
into account the singularities in the lower half plane, we introduce auxiliary functions 

i m (p),f!' 



fii(p)i fiP(p) satisfying the following two conditions: 



- / dpe-* z f^(p) = z3e 



Li 



j -2-Kliz 



dpe~ pz fty(p) =iz J e 
fiP(p)i fij\p) is pure imaginary when p is real. 
Especially, for non- negative j, fif\p) fu\p) are given by 



fli (P) = 7T- 



1 



2m \[p- 2m)i +1 + (p + 2my +1 

j\ ( i (-iy +i 



2m \(p- 2niy +1 (p + 2ixiy +1 



And fi 1 ^ is 



fi,Up) = ~ 7T Hp ~ 2m ) ~ ln (P + 2 ™)] 
In terms of these functions Vq is given by 



^/ff(l>) + 7(?§-CA)/&M 



n=0 



2n + 1 



(26) 
(27) 

(28) 



r) \ 2n+l 

■""'If) (29) 
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where the constant M depends on the choice of the Riemann surface of the logarithm. 
In a similar way, one finds 



± U (p) = \ (c B f[[\p) + ca^Up) + ffSip) 

i 00 / 
= M ' + 4^E (cB(2n + 2 ; 



, ir fP_\ 2n+1 



n=0 

00 _ 



2n + 1/ v ' V2tt/ 

n=0 

where M' depends on the choice of the Riemann surface. 

In order to check the validity of analytical guess (J29)) and (|3U|) about singularities 
and to evaluate ca, cb, we derive power series solutions of (j2"Uj) numerically and compare 
them with the power series expansions of the R.H.S. of (J29*j) an d (|3*U|) . By substituting 
the power series expansions: 

X V j 

^o(p) = ^ a " p " ^b(p) = XXp 71 ; a ° = °' &0 = 

n=0 n=0 

into (J2"U|) and comparing term by term, one can determine the coefficients a n , b n 
recursively. They have the following asymptotic forms. 

417 

-m 2 n+i(-l) n (27r) 2n+1 -> 0.27893(2n + 2) - — -, as n^oo 

{2n + 1) 

M 

= A l (2n + 2) - 



(2n + 1) 

417 

±*2n+i(-l) n (2vr) 2n+1 -> - 0.27893(2n + 2) + — -, as n^oo 

{2n + 1) 

-ia 2n =0, Vra 

±^ 2n (-l) n (27r) 2n -> 0.87628 = A 3 (31) 

The numerical evaluations of At, A3 are quite robust. On the other hand, A2 is sensitive 
to an error of A\ since the coefficients of 1/ (2n+ 1) are fitted after subtracting the leading 
terms. However, we observe that the coefficients of l/(2n+l) in a 2n +i{— l) n (27r) 2n+1 and 
^2n+i( — l) n (27r) 2n+1 are always the same. This observation suggests that the coefficients 
of 1/z in $10, \I/io are identical, which implies 

A2 = 2A = ^ = ^ (32 ) 
m 24m 6 v ; 

The present numerical estimation gives a consistent value of 0.417 with this relation. 

According to (|2H|) . (J3U|) and (pTTjl . the asymptotic expansions of a2 n +i(— l) n (27r) 2ra+1 

and &2n+i(~ l) n (27r) 2n+1 with respect to n should have the same leading terms. This 

is indeed the case. Moreover, A3/A1 should be equal to ir and we have an excellent 

agreement: A3/A1 = 3.14158- • •. Therefore, we can conclude that our guess about the 

poles and branch cuts in the p-domain is valid. Then the constants ca, cb are given by 

ca = iirA2 and cb = i47r 3 Ai. 
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So far, we have considered the first three dominant terms in Vq, Uq. 
From the asymptotic expansions of (|25jh one finds that Vq, Uq involve /y (p) 
(j < -2) and /^(p) (Z > 2). As easily seen, one has the following rough estimations: 

OO OO j 

fl>) ~ E ^^f" w < -2) . /,» ~ E p^f" c* a -i) 

n=0 v ; n=0 v ; 

Therefore, if [/y ] n denotes the absolute value of the coefficient of p n , one has the 
following relation for large n: 

[/ffin » [/ff ]„ » [/ffjn » [/£>]„ " " " • (33) 

As a result, neglected terms are expected to be sufficiently small in the present numerical 
estimation. 

In short, in this subsection, we have constructed the Borel transforms of <&q, 
following the procedure discussed in the previous subsection. The result is summarized 

as 

V (p) =^A 1 ff 1 \p) + ^A 2 f[%{p) 

U (p) = ± (l^AJ^ip) + 2tt 3 A l fg\p) + nA^Up)) (34) 
3.3. Borel Transforms of $oi an d ^oi 

Let <&qj(z), ^Qj(z) be the sum of negative powers of z in the asymptotic expansion of 
Q>oj(z), ^oj(z) and Vj(p), Uj(p) be the Borel transform of v&j^ in the sector including 
z = — oo, then Vj(p) provides the following terms in another sector including z = +oo 

„ oo 

- / dp e-*>* Viip) = J2 ®n 3 (z)e- 2mnz (35) 

T n=l 

In addition, we define $ n (z) = S^=o ^^WO 2 ) f° r later use. Similar quantities \& n j and 
\P n are introduced for 

As the construction of Vq, ?7q from $io, ^xoj Vi, Z7i can be obtained from the 
solutions $n, of the following equations: 

z z 
A* u (s - 1) = ± ^^^y^^e^W ± ^io(^oi(^) e ^ (,) (36) 

Indeed, from this equation, asymptotic expansions of $n, in z are obtained and, 
by exactly the same procedure as in the previous subsection, Vi, U\ are determined as 
(for more detail, see |Appcndix B| ) 

V l (jp) = ± (%^B 2 {k - l)fg\p) + An 3 (B 3 (k + 1) =F B.ktc) f[? {p)) 

Ul{jp) = _ 8 n 4 B 2 (k - l)f$(p) + An 3 (B 3 (k + 1) =F B.kQf^ip) 

-4 7 r 3 J B 4 (fc-l)/i ( ; ) (p) (37) 
And B 1 , B 2 , B 3 and B A can be numerically evaluated as A 1 , A 2 and A 3 : 

B l = 0.7303, B 2 = 0.007399, B 3 = 0.03651, B A = 0.04649 . (38) 
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In this subsection, we study the general structure of the inner solution 
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Yl'jLo^^iji 2 ) ^i( z ) = Yl'jLo a ^iji. z )i which satisfies the following equation: 



A$i(*) = -<7*i(-2;)cos(*o(«)+Moo(*c + ^)) 
A^x(z - 1) = ka$i(z) cos (® (z) + v 00 {t c + az)) 

For this purpose, it is convenient to introduce two solutions of the above 
admitting the expansions: 



(39) 
equation 



oo 



(■A)/ 

y 1 



£>(*) 
*£>(*) 



where $^ (2) and $ x - 



z) are uniquely determined by posing 



coefficient of - in $\.'[z) 

z J 



coefficient of z in ( : ; 



. 



2 



(40) 



Then, because of the linearity of (J39|) . one has the following expressions 



$10 



(B) 



A $S^+A $ 10 



10 



j=0 j=0 



(41) 



which leads to 



$1(2) = (A + ( xA 1 + -- 



+ (A + o-Ai + • • •) ($1? + } + • ■ ■) . (42) 

00 

The coefficient A = o~ n A n is referred to as the Stokes multiplier in ^U] • Let 

n=0 

00 

A = <r n A n , then one finally has 



n=0 



$1 



A 



(A) 



10 



11 



+ A 



(B) 



10 

+ (TV 



11 

(B) 
11 



The first two terms of the Stokes multiplier can be estimated as follows: 
dpe' pz {V (p) + aV 1 (p)} 

dpe' pz {U (p) + aU 1 (p)} 



A 



+ A 



(A), 



-2-7riz 



±A z + a(±A x z ± ^A 2 2 ) 
A 2 + a (A lZ - ^A (z 2 + z)) 



-2niz 



(43) 
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where the contour integral is evaluated with the aid of ()34|1 and (|37|). This leads to the 
following relations 

A = i4tt 3 A 1 = iA%it A B 2 = i24rt 3 B 4 

A 1 = 4ir 3 (B 3 (k + l)TB 1 kt c ) (44) 

Eq. (f4*4"|) implies two linear relations among Ai, B 2 and B4, which are satisfied by the 
present numerical estimations rather well: 

M „ , A 1 



12nB 2 



0.99998 ~ 1 



6ttBa 



0.99997 ~ 1 . 



4. Matching of Inner and Outer Solutions 



4-1. Matching at a Singular Point 

In this section, solutions of the outer equations are constructed and are matched 
with the inner solutions. We first consider the contribution from the singularity 
t\. Corresponding to the expansion of the analytically continued inner solutions: 



Z^n=0 ' * 



Y^=o 2 ™ 2 , the original solutions v and u acquire new 



j n =0 ^ n 

terms in a sector Re t > Re t\ of the ^-neighborhood: 

v{t) — v (t, a) +Vi{t, a)e * + v 2 {t,a)e ° 
u(t) = u {t, a) + ui(t, o)e~'^ Lt + u 2 (t, o)e~^ 



Because of e 2mnz oc e n e~ 



+ ■ 
+ ■ 



(45) 



"* with e = e~°^, v n and u n are of order of e n . By 
substituting (|45j) into (Q) and comparing term by term, we obtain the equations for v\ 
and u\\ 

V\(t + 0) — V\(t) = — cru i{t) cos Uo(t) 

Ui(t + a) — U\{t) = ko~Vi(t + a) cos v (t + a) (46) 
Its solution is uniquely determined by the matching condition: 



vx(t,a) 
Ui(t,a) 



-2-Ktz 



t=t\-\-az 



A 



+ A 



10 

10 

^\z)+a^ ) (z) + 



*)+• 
z) + - 



z)+- 



or equivalently 



v\{t\ + <72, a) 



A 



10 



(A). 
11 1 



z) + 



+ A 



10 
(B) 
10 1 



2) + ^? {z) + 



(48) 



The matching condition suggests the following expansions: 

00 00 
vi(t,a) = o->e—J2vm(t)o- n , u x (t,a) = ^e~ £)u ln (*)<7 n (49) 



n=0 



n=0 
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where j is an integer. Then equations for t> 10 , u w ■ ■ ■ are given by 

u io(*) = - u w(t) cosuoo(t) 
u' w (t) = kv w (t) cos v 00 (t) 

v'nit) + 2 U ioW = cos u oo (t) + u w {t)u i{t) sin m 00 (t) 

«'uW - 2 U io(*) = fe ( u n(*) cos w 00 (t) - Uio(*)«oi(*) sin u 00 (0) 
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(50) 

(51) 
of 



First, we solve (jfiOJ) . Since it admits two linearly independent solutions: 
and the following 

L ^2 J 

- rj- [ sinh Vkt cosh V%f + Vkt ) - k 3/2 tanh Vkt 



x 2 (t) 

V2(t) 



x 2 

-xi(t) 
4fc(l -A;) [ 2^ 



4fc(l - jfe) 



2^ 



f sinh cosh \/A;£ — Vkt] H coth Vkt 



one has 

~v 10 (t) 
u 10 (t) 

with the matching condition 

v w (az + h) 
u 10 (az + ii 



Zl(t) 



+ 6 



x 2 (t) 



(52) 
(53) 



a J < a 



X\{az + 1\) 
y\{az + ti) 



Ar 



An 



+ A 



x 2 (cr^ + ii' 



*S?(*) 



dom. 



(54) 



where the subscript cr° indicates to take the terms of order a and dom. stands for the 
largest part for z — > oo. 

Because L.H.S. of (J54|) starts from z, one should have j = —1. Then, the coefficients 
a and b are determined by the requirement that (|5~4*|) admits well defined o — > limit: 

_ ^(fc-l^ + ^ + fe) 2A 
a " A ° 4ifc(Jfe - 1) ' 
Therefore, up to the 0th order in a, we have 

vi(t) 
Ul {t) 



2A 2^£i h(k - l) 2 + (1 + k) 
— — e ct 



i a 



8k(k-V 



yi(t) 



x 2 (t) 



(55) 



Similarly, up to the first order in a, one obtains 
v x {t) 



Ul (t) 



2A^ 2^ [Uk- l) 2 + (l + fc) 
— : — e °- 



?fc(fc 



Xi(t) 

y 1 {t)+ay[{t)/2 
x 2 {t) 

y 2 (t) + ay' 2 (t)/2 



56) 
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where AW = A + aAi is the first order Stokes multiplier. 

In the t 2 - n eighborhood, the unperturbed solutions v , u acquire the following terms 
in a sector Re t > Re t 2 . 



vi(t) 

Ul{t) 



2A<® 2^ (t 2 (k- l) 2 -(1 + fc) 
e ct 



ia 



8k(k 



x x (t) 

y 1 (t) + ay' 1 (t)/2 
x 2 (t) 

y 2 (t)+ay f 2 (t)/2 



where A^ is the first order Stokes multiplier arising from the singularity t 2 . 



57) 



4-2. Contributions from Other Singular Points 

So far, contributions from t\ and t 2 have been considered. Here, we study contributions 
from other singular points. 



Let t (n) denote the singular points satisfying, |Im[t^ r 



^g^. Then terms 



arising from t^ are of order e 2n+1 , where e = e . Since (vo(t), uo(t)) is bounded and 



(v n (t),u n (t)) 



[n > 1) for sufficiently large t, one can neglect the terms arising 



from (n > 1) even for sufficiently large t. 

On the other hand, terms arising from t± may produce new terms when t passes 
through t 2 - n eighborhood, but they are negligible as they are higher order with respect 
to e. Therefore, the overall contribution of the singular points is a simple sum of the 
contributions from t± : t 2 , t\, t\ and one has the following solution of the unstable 
manifold in the whole time domain: 



v u (t) = Voo(t) + cr 2 v 02 (t) 

~4AW 2^ 
e <> 



+ S_(t)Re 
+ S + (t)Re 

Unit) = u 00 (t) 

+ S-(t)Re 



tjjk- 1) 2 + (1 + fc) 
ia \ 8fc(fc-l) 

4A( 2 ) 2« it , ( t 2 {k- l) 2 - (1 + fc) 
ia 6 ° V 8fc(fc- 1) 



xi(t) - x 2 (t) 



Xi(t) - x 2 (t) e 



^ + a 2 u 2(t) + a 3 Qt4(*)-^(*) 

4AW gTTgi r 
— : — e CT < 



ti(fc- 1) 



+ 5+(t)Re 



2(7 



4A( 2 ) 2^ 

e = 

ia 



8fc(fc- 1) 

M *- 1)! - (1 + fc) 



8fc(fc- 1) 



S/ 2 (*) + 



3/2 _*»*' 



(58) 
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where S±(t) = S(t±tf) and S(t) stands for the step function, and AW, A( 2 ) denote the 
first order Stokes multipliers from the analysis of ti, t2, respectively. 

4-3. Comparison with Numerical Calculation 

Here we compare the analytical solution fl^Sj) obtained in the previous subsection with 
the numerical calculation. For this purpose, it is convenient to rewrite the solution as 



v u {t) = v {t) + 2S(t) 

+ 2S(t + 2tf) aia(t) cos 



an(t) cos h On(tj sm 

a a 



2vr(t + 2tf 



6l2(0 sin 



. 27r(t + 2tf ) 



a 



u u {t) = u (t) + 2S(t) 



2nt / \ • 2vrt 

«2i {t) cos h 621(f) sm 

a a 



+ 2S(t + 2t\ 



a 22 (t) cos 



2n(t + 2t[ 



a 



+ 622(f) sin 



2vr(t + 2tf ) 



where the auxiliary functions are given by 

an(t) — (K + K) (ax^t) - x 2 (t)) - K'^x^t) 

K[ (ax^t) - x 2 (tj) +(K + Ki)^i(i) 



b n {t) 



a 12 (t) = (Kt - K) (axi(t) + x 2 (t)) ~ K[0Xi(t) 



b l2 (t) = - K[ (axtit) + x 2 {t)j + (K — K x )fJx x {t) 
a 21 (t) = (K + K x ) [a fat) + |&(*)) - fat) + |&(t) 

-K[(3 fa^ + fatj) 
b 21 (t) = K[ {a fat) + a -y[{t)) ~ fat) + } 

+ {K + K X )P fat) + ^y[(t)) 
a 22 (t) = {K x - K) [a fat) + |#(t)) - fat) + |&(t) 

-ifi/? fat) + ^m) 

- k[ { a fat) + a -m) + fat) + ^m) } 

+ (K-K 1 )P fat) + a -y[{t) 



(59) 



&22W 



(60) 

The functions with tilde are defined as /(f) = f(t + tf) and the constants a, (3, K 1} K[ 
are given by 

87^! 



K 



a 



Kt = -8ir 3 Bxt{e m 
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a 



87r 3 (fi!Hf -B 3 (k + l))e 
tf(k - l) 2 + (l + k) 



Ts/k 



t[(k-l) 



(61) 



8Jfe(Jfe-l) ' r 8k 

where tf and t[ are the real and imaginary parts of t\, respectively. Eq. (|59jl indicates 
that, when t exceeds t — 0, exponentially growing oscillatory term appears and, when 
t exceeds t = — 2tf(> 0), another exponentially growing term is added. These terms 
should describe heteroclinic tangles. 

This is indeed the case. Fig. HJ Fig. [S] shows an analytical solution of the unstable 
manifold from (tt,0). One can see a large-amplitude oscillation near (— tt, 0). As shown 
in Fig. |H Fig.0 the analytical solution (JoTJj) (solid curve) well reproduces the oscillation 
in the numerically calculated orbit (dots). The dots show the time evolution of an 
ensemble consisting of 500 points. One can see the stretching of the ensemble is also 
well reproduced by the approximate unstable manifold. When only the leading order 
terms of the inner equation are retained, the agreement between the analytical and 
numerical results is not so good (cf. Fig. EJ) and, thus, the terms in the inner solution 
of order a are important. 

Before closing this section, we give an approximate stable manifold (v s (t),u s (t)). By 
the similar approach to construct the unstable manifold, we have the similar expressions: 



v s {t) = v (t) - 2S(-t) 

- 2S(-t + 2tf 
u s {t) = u (t) - 2S{-t) 

- 2S{-t + 2tf 



2-nt , . . . 2-nt 

an (t) cos h On yt) sin 

a a 



a 12(2) cos 



2vr(t + 2t[ 



a 



+ b 12 (t) sin 



2n(t + 2t[ 



a 



2nt , . . . 2nt 

a 2 i{t) cos h b 2 i{t) sin 

a a 



a 22 (t) cos 



2n{t + 2tf ) 



a 



+ b 22 (t) sin 



2ir(t + 2t? 



a 



(62) 



5. Reconnection of Unstable Manifold 

We remind that the separatrix of the continuous-limit equation is the unperturbed 
stable/unstable manifold. As mentioned in Introduction (cf.Fig. [TJ, the separatrix 
changes its topology depending on the parameter k (reconnection). When < k < 1, 
there appears a separatrix connecting (tt, 0) and (— n, 0). As k — *■ 1, it approaches 
the union of two segments: the one connecting (tt, 0) with (0, 7r) and the other (0, tt) 
with (— 7r, 0). Both segments are separatrices at the parameter value k — 1. The 
corresponding topological change occurs for stable/unstable manifolds, which will be 
discussed in this section. As shown in the previous section, since the approximate 
stable and unstable manifolds are related with each other by simple symmetry, it is 
enough to discuss the topological change of the unstable manifold. 
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[v - ?r) x 10 fi 



Figure 4. The analytically constructed unstable manifold (solid line) near (— it, 0) and 
time evolution of the ensemble (dots). Inset shows the overall view of the analytically 
construced unstbale manifold, a — 0.35, k = 0.85 




Figure 5. The analytically constructed unstable manifold (solid line) near (— w, 0) 
and time evolution of the ensemble (dots), a — 0.35, k = 0.85 

Fig. [7[ Fig. |H1 Fig. El Fig. El show the unstable manifolds when a = 0.35 and 
k = 0.2, 0.5, 0.85, 1.0 — 10~ 9 , respectively. The overall structure of the unstable 
manifold except near the fixed point (— ir, 0) is very close to that of the separatrix. 
Near (— 7T, 0), the unstable manifold acquires an oscillatory portion and the slope of this 
portion becomes steeper as k — > 1. The behavior of the slope can be understood from 
the asymptotic ratio between the additional terms: 

lim — - — = —y/k cr . 

i^+oo Vi(t) + 7T 2 

It is interesting to see that, when k is very close to unity (i.e., k = 1.0 — 10 -9 ), there 
appears an additional oscillatory portion in the unstable manifold near (0, n) (cf. Fig.[SJ). 
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Figure 6. The analytically constructed unstable manifold (solid line) near (— w, 0) 
and time evolution of the ensemble (dots) . Only the leading term of the inner equation 
is taken into consideration, k — 0.85, a — 0.35 



3 



u 



-3 




-3 3 

v 



Figure 7. The analytically constructed unstable manifold, a = 0.35, k = 0.2 

As mentioned before, at k — 1, two segments from (tt, 0) to (0, it) and from (0, 7r) to 
(— 7r, 0) are separatrices of the continuous- limit equation. Thus, the perturbed unstable 
manifold at k = 1 starting from (jr, 0) should have an oscillatory portion near (0, it). 
The oscillation near (0,7r) of the unstable manifold shown in Fig. can be considered 
as the precursor of the oscillation in the unstable manifold at k = 1. 

The origin of these behaviors are summarized as follows. When k is not very close 
to 1, splitting term is exponentially small compared to the unperturbed solution and 
can be negligible for not too large t . But this term becomes dominant for sufficiently 
large t because 22, 2/2 — ► 00 as t — » 00. Thus, a large-amplitude oscillation appears 
near (— 7T, 0). This is not the case when k is very close to 1 because exponentially small 
terms with respect to a are proportional to t^t. Therefore, the additional terms are 
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Figure 10. The analytically constructed unstable manifold, a = 0.35, k = 1.0— 10 
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v x If! 1 

Figure 11. The analytically constructed unstable manifold and the oscillation near 
(0,tt). a = 0.35, k = 1.0- 10" 9 

not too small compared to the unperturbed solution near (0, it) and we can observe the 
oscillation near (0,7r). We further remark that two oscillatory portions in the unstable 
manifold for k ~ 1 come from two singularities ti and t 2 in the complex time plane. In 
other words, the two sequences of singularities are necessary to exhibit the reconnection 
of stable/unstable manifolds. 

6. Conclusions 

With the aid of ABAO (the asymptotics beyond all orders) method, we have derived 
analytical approximations of the stable/unstable manifolds, which agree rather well 
with those obtained by the numerical iteration. When k is not close to 1, the perturbed 
unstable manifold starting from (rr, 0) exhibits highly oscillatory behavior near the fixed 
point (—7T, 0) and its overall structure far from (— tt, 0) is almost the same as that of the 
unperturbed unstable manifold. However, when k is very close to 1, we can observe a 
oscillation near (0, tt), which is considered to be the precursor of the heteroclinic tangle 
in the unstable manifold at k = 1. In this way, even when the heteroclinic tangle exists, 
the unstable manifold smoothly changes its topology as the change of the parameter k. 

Contrary to the systems studied so far by ABAO method, the unperturbed solution 
of the Harper map have two sequences of singular points in the complex time plane and 
the interference of the contributions from them might be expected. In this paper, the 
approximate stable/unstable manifolds are constructed just adding the contributions 
from the two sequences, but they agree quite well with the numerically constructed 
unstable manifold. In a higher order approximation, we need to take into account this 
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interference. This aspect will be discussed elsewhere. 
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Appendix A. Case of k > 1 



In this section, we consider the case of k > 1. Put 



t = kt, a = ka, u(t) = v 



then v, u are the solution of the Harper map with parameter 1/k: 



v(t + cr) — v(t) 
u(t + a) — u(t) 



asmu(T\ 



a 



sinv (t + a) 



(a.i; 



(A.2) 



Appendix B. Borel-transformed first order solution of inner equation 

The aim of this appendix is to calculate the Borel transforms of the solutions $q 1; \E' 01 . 
In order to investigate the behavior near t — t\ and t = t 2 in a unified way, we introduce 
$oi 0) and * i(*) by 

(B.l) 



'oi = $oi > ^oi = *oi ' 
The Borel transforms V\ , U\ of $oi , ^ oi satisfy 

(e- p - 1)U X 



-— V(p) + \) + ^j^(g(p) - 1) + Ut * g{ P ) 



;i - e p )f/i = ^Af(p) + Vi * f(p) 



where 



B 



2 

z 



f(p), B 



3 =F«*o 



g(p) 



The power series expansions of V\, U\is defined by 



Vi 
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n=l 



n=l 



(B.2) 



(B.3) 
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c n = ±kt e $» + + , d n = ±ktd$ + kd^ + 43) (R4) 

where ch\ dn are independent of k. By substituting (|B.4|) into ()B.2|) . we numerically 
get the following estimation. c n , d n as follows: 

A2n + 2)(2n + 1) , , „ u 2n + 1 

c 2n = B 2 (k - l) y - J n ) 1 . ; - (±B 1 kt c - B 3 (k + 1))- 



(27r) 2 «(-l)« v 16 JV // (2tt) 2 ™(-1)' 

~ , 1 ,(2n + 2)(2n + l) , , „ NN 2n + 1 

d 2n = B 2 (l - k) - ') . ' - (±B 1 kt c - B 3 (k + 1)) 



d2n+l = ~ B A (k - 1 



(27r) 2 ™(-l)« v c jv // (2tt) 2 ™(-1)" 
(2n + 2) 



(27r) 2 " +1 (-l) n 
This estimation gives (pTTjl . 
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